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Abstract 

We use Lie point symmetries of the 2+ 1-dimensional cubic Schrodinger equation 
to obtain new analytic solutions in a systematic manner. We present an analysis 
of the reduced ODEs, and in particular show that although the original equation is 
not integrable they typically can belong to the class of Painleve type equations. 

1 Introduction 

The purpose of this article is to find new explicit group invariant solutions of the cubic 
Schrodinger equation (CSE) 

iV> t + AV> = aoMV (1-1) 

where ip(t,x,y) is a complex-valued function of its arguments, ao is a real constant and 
A is the Laplace operator in M 2 . 

Eq. (I1.1J1 plays an important role in many branches of physics. For instance, in nonlin- 
ear optics it describes the interaction of electromagnetic waves with different polarizations 
propagating in nonlinear media such as an isotropic plasma. Another instance where Eq. 
fll.ljl arises is water waves in a deep fluid. 

In one space dimension, the CSE is known to be integrable via inverse scattering 
transformation. In two space dimensions this approach does not work. Therefore, our 
approach will rely on the use of powerful methods of symmetries of differential equations. 
The Lie group G of local point transformations leaving the equation invariant is well- 
known in literature. We shall use subgroups of G to reduce the equation to a system 
of algebraic or (second order) ordinary differential equations. At this stage, it is of vital 
importance to perform a classification of subalgebras. The reduced equations will be 
compatible with specific types of boundary conditions, such as translational invariance 
or a cylindrical symmetry. Then, they will be solved whenever possible, in particular 
when they have the Painleve property. This approach is basically similar to that of a 
generalized non-linear Schrodinger equation involving cubic, quintic equations as special 
cases in three space dimensions (n = 3) studied in a series of papers [H 121 E] . For a vector 
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version of (jl.lj) . a similar analysis for two space dimensions (n = 2) and two or three wave 
functions (N = 2,3) is given in Ref. jl]. The present article will complete the results of 
[4J to include the scalar case, namely one wave (N — 1). 

We organize the article as follows. In Section 2, we review the symmetry group G, its 
algebra L and subalgebras of the Schrodinger equation. In Section 3, we utilize subalgebras 
to perform reductions to second order ordinary differential equations (ODEs). In Section 
4 we construct solutions of the reduced equations whenever possible and thereby invariant 
solutions of the original equation. In particular, we show that for specific values of group 
parameters, we can pick out equations having Painleve property among generically non- 
Painleve ones. Finally, in Section 5 we make some concluding remarks. 



2 The symmetry group and its Lie algebra 

By a symmetry group G we shall refer to the local Lie group of point transformations 

i=T g (t,x,y,ip,ip*), x = X g (t,x,y,^,ijj*) } y = Y g (t,x,y,ip,ip*), 
■tp = ^g(t,x,y,ijj,i/j*) 

such that i/;(t,x,y) is a solution whenever ip(t,x,y) is one. Here star denotes complex 
conjugation and g is the group element. 

The computation of the symmetry group is quite straightforward (for details see for 
example |5|). To this end, we look for the infinitesimal symmetries that close to form a 
Lie algebra under commutation. The Lie algebra L of the symmetry group G is realized 
by vector fields of the form 

v = rd t + £d x + 7]d y + cpd^ + ip*d^,* . (2.2) 

The coefficients r, £, rj, <p, (p* are functions of t,x,y,ip,ip* to be determined from the 
invariance condition 

pr (2) v(£) = 0, (2.3) 

E=0 

where E is the equation under study and pr^ is the second prolongation of the vector 
field v. This condition provides us a system of determining equations for the coefficients 
of v. Solving this system we end up with a basis (depending on integration constants) of 
vector fields for the symmetry algebra L. We exponentiate them, namely calculate exp L 
to find the symmetry group G. 

An inverse problem of determining the most general form of Schrodinger type equations 

if) t + F(x, y, t, ^, Vi, Aj, rl>*j) = 0, i,j e {x, y} (2.4) 

invariant under the Schrodinger and its subgroups was attacked in . An investigation of 
this problem was reduced to that of classifying inequivalent realizations of the Schrodinger 
algebra. In particular, the CSE appears to be among those equations allowing the same 
Schrodinger algebra set) (2) as its symmetry algebra. Hence, we simply resort to the 
results of |J] for recovering a suitable basis of the algebra written in phase and modulus 
coordinates. In passing, we mention that symmetry breaking interactions for the time 
dependent Schrodinger equation were studied in [7j for n = 1 and in [S] for n = 2. 
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We introduce the moduli and phases of ip by setting 

ip = Re**, < R < oo, < <p < 2ir. (2.5) 

and rewrite our equation (jl.lj) in terms of R and as a coupled system of real equations 

-R(f>t + AR-R\V<p\ 2 = a R 3 , (2.6a) 
Rt + RA(p + 2VP ■ V0 = 0, (2.6b) 

where R and <fi are functions of t, x, y. 

We summarize some facts on the structure of the symmetry algebra and its important 
subalgebras. 

• The nine-dimensional Schrodinger algebra set) (2) has a Levi decomposition 

scfj(2) = fj 2 e{sl(2,R) ©so(2)} 

~ {P x , P 2 , Pi, P 2 , E} e{T, C, P>, J} (2.7) 

with nonzero commutation relations 



[Pi, Pi 


] = 


[P 2 ,P 2 ] 


= E/2, 




= -Pi, 


[J, P 2] 


= -Pi, 


[J,B 1 ] 


= P 2 , 


[J, Pi] 


= P 2 , 


[T, B 3 ] 




[D, Bj] 


= Bj, J = 1,2 


[D,Pj] 


= 


[C, Pj] 


= -B 3 , j = 1,2 


[T, D] = 


= 2T, 


[T,C] 


= D, [D, C] = 2C. 
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Here <Ef- denotes a semi-direct sum and f) 2 is a Heisenberg algebra. 

The subalgebra {T, P x , P 2 , B x , B 2 , J, E} generates the extended Galilei group (trans- 
lations, proper Galilei transformations and constant change of phase), D and C 
generate dilations and nonrelativistic conformal transformations. 

A convenient basis for the standard Schrodinger algebra is given by the vector fields 

T = d t , P X = 8 X , P 2 = d y , E = d 4> , (2.9a) 
B x = td x + x/2d^, B 2 = td y + y/2d 4> , J=-yd x + xd y , (2.9b) 
D = 2td t + xd x + yd y - Rd R , (2.9c) 

C = t 2 d t + xtd x + ytdy - tRd R + ^{x 2 + y 2 ^. (2.9d) 

The transformations corresponding to the Galilei-similitude algebra are global and 
given by 

t = e x (t-t ), 

x = e A/2 {cosa[(a; - x ) + v x (t - t )) + sina[(y - y ) + v 2 (t - t Q ))}, 

y = e A/2 {-sina[(x - x ) + v x (t - t )} + cosa[(y - y Q ) + v 2 (t - t )}}, (2 ' 101 

R = e - x ' 2 R, 4> = <f)+ -{v x {x - x ) + v 2 {y - y )} + ^{vf + v 2 2 )(t - t ) + 5, 

where t , Xq, y$, v x ,v 2 , a, X, S are group parameters. 
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The conformal transformations generated by C are local and given by 

t x _ y 
t = ~i 1, x = - -, y 



1 — pt 1 — pt " 1 — pt 

p(x 2 + y 2 ) 



R = (l-pt)R, 



r 4(1 -pt) ' 
provided 1 — pt ^ 0. Above, p is group parameter. 

The transformations (j2.1()|) and (j2.11j) can be composed to obtain the transforma- 
tions of the Schrodinger group. 

The CSE is also invariant under time reversal and reflection of the space coordinates: 



T d : t-* -t, x -> x, y -> y, ip ^ip*, 

R x : £->•*, x -> -x, y-^y, i>-+$, (2.12) 

R y : t — > t, x — > x, y — > —y, if) — > if). 

The elements T, R x , R y generate a discrete finite group £?£>. 

Comment: The conformal transformations generated by C do not generalize to higher 
dimensional cubic Schrodinger equations. We mention that in a higher dimensional space 
n, the nonlinear Schrodinger equation with a critical exponent 4/n |Hj 

iip t + Aip = \ip\ 4/n tp, xeW 1 

is conformally invariant. In this case, the conformal generator C and the corresponding 
local group action on (x, t, R, (f>) are given by 

n + I 1 2 n 

C = t 2 d t +tY^Xjd x . - ^Rd R + ^«9 , \x\ 2 = (2-13) 

3=1 3=1 



X J 



1-pt' i3 ~\-pt' ^" 1 '---' n 

1 ' (2.14) 



R= (l-pt) n / 2 R, 



p\x\ 



4(1 -pt)' 



provided 1 — pt ^ 0. The induced action of C on a solution if> (x,t) produces the trans- 
formation formula 



iP(x,t) = (l+pt)-"/ 2 expf * m - V/— ^— ,- 

V y> ^Uil+p^J^Kl+pt'l+pt 

which means that ip(x,t) is a solution whenever ipo(x,t) is one. 



; ' xeR n 
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2.1 Subalgebras of the Schrodinger algebra scf)(2) 

To perform symmetry reduction for Eq. (jl.lj) we need a classification of low-dimensional 
subalgebras of the extended Schrodinger algebra L = scf)(2) into conjugacy classes under 
the action of the group of transformations G leaving equation (jl.lj) invariant. The group 
G has the structure 

G = Go (X Go, 

where (X denotes a semi-direct product. The invariant subgroup G is the connected com- 
ponent of G, the Lie algebra of which is (j2.7j) and Gp is the discrete component defined 
by (j2.12j) . In order to reduce Eq. (jl.lj) to an algebraic equation, an ordinary differential 
equation (ODE) or a partial differential equation (PDE) in two variables we need to know 
all subalgebras Lq C L of dimension dimLo = 3,2 and 1, respectively. Moreover, the 
corresponding subalgebras should not contain the element E alone as long as we are in- 
terested in invariant solutions. Otherwise, the subgroup invariants would be independent 
of the phase <fi and this in turn would not define locally invertible transformations to ip 
or ip*. 

A complete list of representatives of all subalgebra classes of scf)(2) are given in Ref. 

ma. 

3 Symmetry Reductions 

A group-invariant solution is one that is transformed to itself by group transformations. To 
find such a solution we choose a suitable canonical subgroup and require that the solution 
be left invariant. In this section we are interested in looking at solutions invariant under 
two-dimensional subalgebras since they lead to reductions to real systems of ODEs. We 
go through each representative subalgebra and find the corresponding reduced system 
and solve them whenever possible, and thereby a solution of the original PDE. One can 
transform these solutions by group transformations to obtain all other possible invariant 
solutions. Note that the ideas we employ here are described more fully in [H 121 IH] • 

In table [T] we list representatives of two-dimensional subalgebras that provide reduc- 
tions to ODEs. We note that we skipped those trivially acting on the coordinate space 
(t,x,y). 

Subgroups that have generic orbits of codimension one in the coordinate space (t, x, y) 
and of codimension three in the total space (t, x, y, R, (p) will lead to reductions to ODEs. 

The procedure for finding invariant solutions is to require that the solution ip(t,x,y) 
be invariant under each subalgebra. This invariance imposes constraints on the solution 
that are expressed by first order linear PDEs. Their solutions imply that the invariant 
solutions should have the form 

V(f , x, y) = f(t, x, y)M(0 exp [i(g(t, x, y) + $(0)], M(£), $(0 e R, (3.1) 

where f(t,x,y), g(t,x,y) and £(t,x,y) (the subgroup invariant) are explicitly known for 
each subalgebra. Below we give the expressions for £ and the wave if>(t, x, y) together with 
the reduced equations for the functions M(£) and in a list. Note that for all of the 
abelian subalgebras, one of the coupled reduced equations, namely reduction of (j2.6bj) 
is directly integrable once. Using this we eliminate to obtain a second order real 
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Table 1: Two-dimensional subalgebras leading to ODEs (e = ±1) 



No 


Type 


Basis 


£2,1 




J + aE,C + T + bE 


L*2,2 




J + aE,T 


-^2,3 




J + aE,T + eE 


-^2,4 




J + aE,D + bE 


£2,5 




T,P 1 


-^2,6 




T + eE, Pi 


L2J 




T + B 1 ,P 2 


-^2,8 




Pi,P 2 


-^2,9 




B U P 2 


^2,10 




J — e(C + T) + a£, P x + eP 2 


-^2,11 




D + aJ + bE,T 


-^2,12 




D + a£, Pi 



ODE for M(£). The non-abelian cases are exceptionally difficult to handle. In the last 
case (1/2,12), w e raise the order of the ODE to make a decoupling possible and hence have 
to deal with a third order ODE. The subalgebra L 2 ,n is the only case that has led to a 
reduction which we have not been able to succeed in treating analytically in any way. 



The subalgebra L 2 1 



i(a8+b arctan t+t£, /4+*(£)) 

£ = --, r-75- , ^ = arctan-, 
1 + ^ x 



The subalgebra L 2)2 



(3-2) 



M 2 $£ = C , C = const. (3.3a) 
^4 + ^)M-f + fM-3 + |M3. (3.3b) 



^(f,x,j/) = M(p)e^ + *^\ 

p=( a ; 2 + 2/ 2 ) 1 /2 ? = arc t a n^, (3 ' 4) 



M 2 $p = C , C = const. (3.5a) 

n 2 r 2 M 

M -,M ■ — ^M~ 3 - — + a M 3 . (3.5b) 
p 2 p 2 p 
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The subalgebra L 2 3 

ijj(t,x,y) = M(p)e l{ae+£hlt+ ^ p) \ 
p = ( x 2 + y 2 ) 1/2 , 9 = arctan -, 



x 



M 2 §p = C , C = const. 

n 2 C 2 M 

M = (e + ^)M +^M" 3 - — + a M 3 
p 2 p 2 p 



The subalgebra L 



2,4 



£ = — 9 = arctan — , 
t x 



M 2 (8$ - l)f = C , Co = const. 



The subalgebra L 



2,5 



i/>{t,x,y) = M(y)e i *M, 



The subalgebra L 2 ,6 



M 2 $ = C , C = const. 

-'0- 



M = C 2 M" 3 + a M 3 . 



il>(t,x,y) = M(y)e i{£t+ ^ v)) , 



The subalgebra L 



2,7 



M 2 $ = C , C = const. 
M = eM + C 2 M" 3 + a M 3 . 



^(t,x,y) = M(0e i ^-V+»(«), 
t = X ~ 2' 



M 2 $ = C , C = const. 

M = -M + C 2 M~ 3 + a M 3 . 
2 



• The subalgebra L 2 8 

ij(t,x,y) = M(t)e^\ (3.16) 



$ = -aClt + d, C , Ci = const. (3.17a) 

M = C . (3.17b) 

• The subalgebra L 2 ,9 

^(t,x,y) = M(t)e l ^+*W), (3.18) 



The subalgebra L 



2,10 



M = C t~2, (3.19a) 
$ = -2aC ** + Ci, C , Ci = const. (3.19b) 



<f>(t,x,y) = — - aearctant + - ^ - ^ + $(0, ( 3 - 20 ) 

= ex-yt 
* 1 + t 2 ' 



M 2 $ = C , C = const. (3.21a) 

M = (£ 2 — ae)M + C^M" 3 + a M 3 . (3.21b) 

• The subalgebra L 2 ,n 

^t,x,y) = ^e<^), 

£ = pexp(-#/a), (3-22) 
p = (x 2 + y 2 ) 1/2 , # = arctan-, 

(1 + a 2 )£M$ + 2(1 + a 2 )£M$ + (1 - a 2 )M$ - 26M = 0, 

(a 2 - 6 2 )M + (1 - a 2 )£M + (1 + a 2 )£ 2 M + 2&£M$ - (1 + a 2 )£ 2 M$ 2 ( 3 - 23 ) 
= a 2 a M 3 . 
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The subalgebra L 2 ,i2 



^x,y) = ^e'(«^) : 



(3.24) 



6 = Co + f x = M 2 , Co = const. (3.25a) 

A, 

2XX - X 2 + - 2a)x 2 - ^X 2 - 4a X 3 -^X-| = 0. (3.25b) 

4 Discussion of the solutions of reduced ODEs 

In Section 3, we showed that for the abelian subalgebras of type 2A\ the reduced pair of 
equations can always be decoupled by integrating the second one once. For the nonabelian 
case of type A 2 , decoupling is possible only for one subalgebra at the price of introducing 
third order derivative. Luckily, it happens to have the Painleve property and its solution 
is expressible in terms of the fourth Painleve transcendent Pjy. 

For each reduction corresponding to subalgebras L 2l i, i = 1, . . . , 10, M(£) satisfies a 
second order ordinary differential equation. All of them can be written in a unified form 
(see Table EJ) 

M = Si(0M + Ai(0M + A 3 (£)M 3 + B Q {i)M-\ (4.1) 
We change the dependent variable to H by putting 

M = VH, H > (4.2) 

so that the Eq. is transformed to 

H = — H 2 + + 2A!!! + 2A 3 H 2 + —. (4.3) 
2H H K J 

To simplify ()4.3|) we apply the transformation 

H = \{£)W(7i), V = V{£) (4-4) 



and obtain 




-,2 \ i ^2 ^ fj2 2\ 2 rf Xrj 2 



w 2 + ^w-\ 



(4.5) 



ff \ 2 ff 



Eq. ()4.5|) can be transformed to one of the canonical equations classified by Painleve 
and Gambier (see jTTj for a list of canonical equations). We shall see that only for special 
values of coefficients they can be integrated in terms of solutions of linear equations, 
elliptic functions or Painleve transcendents. Otherwise, they would not have Painleve 
property which means that they can not be integrated. 
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Table 2: Coefficients of (jUJ) 



No 




A 3 


5 


B 1 


-^2,1 
-^2,2 
-^2,3 
-^2,4 
-^2,5 
-^2,6 
-^2,7 
-^2,10 


6/(40 + a 2 /(4£ 2 ) + 1/16 

«7p 2 

5 _|_ a 2 / p 2 
6/(80 + a 2 /(4a- V64 


£ 

CP 

e - ae 


ao/(40 
a 
a 

ao/(40 
a 
a 
a 
ao 


CUP 2 
CUP 2 

cum 2 ) 

ci 
ci 
ci 
ci 


-ve 
-1/p 
-1/p 
-m 








4.1 Solutions using the subalgebra Z^i 
Eq. (jOJ for (|3~3b"J) is 



J_ _6_ a 2 1 A A 2 A x 
7^2f + 2^ 2 + 8~A£ + 2A 2 ~A' 



(4-6) 



This equation may be transformed to either PXXXIII or PXXXIV among the canonical 
equations ([HJH2])- F° r completeness, we quote them here: 

PXXXIII: W = ^W 2 + 4W 2 + 7 W-^W" 1 , (4.7) 
PXXXIV: W = ^W 2 + 4 7 W 2 -7?W-iw- 1 . (4.8) 

Comparing the coefficient of W 2 and W in ()4.(ij) with (J4.7)) and ()4.8)1 we see that A and r\ 
must be chosen as 

A = A C 1/3 , V = Voe /3 - (4-9) 
For this choice of A and 77 ()4.6|) becomes 



w = —w 2 + ^w 2 + is^fw- 1 + 

2W 2 V 2 l A V 



(9a 2 -1) o 96 9 4 
'-if H =ri H -,r\ 



W. (4.10) 



Inspecting (14. 7j) and ()4.8j) . we see that the coefficient of W must be a constant or rj. For 
Eq. ()4.1()j) . in the coefficient of W, it is possible to remove the term vf 2 by choosing 
a 2 = 1/9 but we cannot get rid of the term rf. Hence Eq. (J3.3bj) cannot be transformed 
to one of the canonical forms, thus it does not have the Painleve property. 

If we consider Eq. fl3.3b|) with Co = 0, similar arguments apply. We mention that Eq. 
(j3.3b|) admits no point symmetries and no constant solutions. 
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4.2 Solutions using the subalgebra L2 2 
4.2.1 Equation (l33bl . with C ^ 

By a suitable choice of A and 77 this equation can be transformed to a standard form which 
has no first derivative and a constant coefficient of the quadratic term. The transformation 
and the corresponding transformed equation are given by 



The above equation has the Painleve property only when a 2 = 1/9. In this case we 
set 7/0 = 9aoAo/8 and integrate once to obtain 

2C 2 

W 2 = 4{W 3 + CW -5 2 ), C = const., 5 2 = — ^ (4.11) 

a A 

whose solution can be expressed in terms of Weierstrass function or Jacobi elliptic 
functions ^3]. By virtue of (|4.2j) we require XW to be positive. Writing the right-hand 
side of fETTIl) as 

P(W) = A(W 3 + CW - 5 2 ) = 4(W - Wx)(W - W 2 )(W - W 3 ), 

where W\,W 2 ,W 3 are roots of the cubic polynomial P(W), we investigate the solutions 
of ()4.11j) for four different cases: 

(i) P(W) has only one root. This is not possible since P{W) does not include the 
term W 2 . 

(ii) P(W) has a double root, say W\ = W 2 ^ W 3 and write P(W) as 

P{W) = A(W 3 + CW - 5 2 ) = A(W - W{f{W - W 3 ). 

— 1 4 12 

This case is possible for C = — 3.2~5s and ^1,2,3 are given as Wi )2 = —2~35a, W 3 = 
(25)3. Integrating for W > W 3 we reach to the solution given below: 

w 2ci . 2 . 1 _ i,l . 1 2 

M = — (— )5p 3 (__ + sec^rja, r = cip3+c 2 , cf = — (a C ) 3 

3 an o 8 , 

r~ r~ 4 - 12 

3 2 1 

$ = tan _1 (y - tanr) - J -r + $ , V = Mexp[7(±-^ + $)]. 

(iii) P(W) has three distinct roots. In this case, roots ^1,2,3 are determined by the 
system 

W 1 + W 2 + W 3 = 0, WxW 2 + W X W 3 + W 2 W 3 = C, W 1 W 2 W 3 = 6 2 . 



We see that two of the roots have the same sign, while the third is always positive. Their 
solutions depend on ordering of the roots. 

S 2 

(a) Wi < W 2 < < W 3 (6) < W x < W 2 < W 3 , W 2 2 < — . 

vV\ 



11 



Integration for W > W3 gives both for (a) and (b) 

W = W 2 + (W 3 -W 2 )[cn(r,k)}- 2 , r = c lP f + c 2 , k 



4 , ,2 Wa-Wi 



W 3 - Wi 



M = p-^^Wf'^ c\ = 9 -^(W 3 - W 1 ), (4.13) 

o 

* = VW=W 1 .[W + <S '<>' ^ = ^exp[ i( ±i 9 + *)] 

with a , Ao > 0. Integrating on Wi < W < Wi we get 

iy = iy 1 + (iy 2 -iy 1 )sn 2 (r,fc), r = Cl p3 +C2 , fc 2 ■- - 



W3-W1 



M = p~ 1/3 (XqW) 1/2 , c l = ^ { W 3 - Wl ) (4.14) 

o 

* = 7w=m I % + *»• * = M -pi ! (±5 9 + *)i 

with a , Ao < for the case (a) and a , \ > for the case (b). 

(iv) Last we investigate the case where P(W) has two complex roots and a real root. 
Let W2 = p + iq, W$ = p — iq, then p, q and W\ are found from the system of equations 

2p + ^ 1 = 0, f- + 2 P Wi = C, p 2 + g 2 = |!-. 

Since W± > 0, integration on W > W\ provides the following solution for ao, Ao > 

jxr A + W1 + (Wi — A)cn(r, k) 
1 + cn(r, k) 



& f d-T ^ 2 3AaoAo 



(4.15) 



^ = Mexp[i(±-0 + $)], A 2 = (p - I-Vi) 2 + g 2 . 

3 

4.2.2 Equation (|3~5b"|l . with C = 
From (j4.11j) We have 

VT 2 = 4W(C + W^ 2 ), C = const. (4.16) 

Depending on the integration constant C, we obtain the following solutions for (jl.ljl : 

(i) C = 0. In this case the solution is written as 

V> = \{-) l, \p 2 " ~ Po)" 1 exp[i(±±0 + $ )], $0 = Const.. (4.17) 

O CLq O 

(ii) (7 = -p 2 < 0. We write Eq. KT^i as 

W 2 = AW(W-p)(W + p). 
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Integration for W > p and —p < W < gives respectively 

r = ^ i (-) V V 1/3 [cn(c 1 p 2 / 3 + c 2 ,A;)]- 1 exp[z(±^ + $ ( ,)i. „ .A > 

3 ao 3 

' ' = ^(— ) V V 1/3 cn (c lP 2 / 3 + c 2 , fc) exp[z(±^ + %)], a t) . A < 



(4.18) 



with c 2 = 9aoAo( 4 " C)1/2 ,fc 2 = 1/2. 

(iii) C = p 2 > 0. The solution for a , A > is given by integrating ()4.16j) 

V- = ^(-) 1/2 p- 1/3 tn(c lP 2 / 3 + c 2 , fc) dn(c lP 2 / 3 + c 2 , fc) exp[z(±^ + $ )] (4.19) 

O Go O 

with c 2 = 9aoA f V2 , fc 2 = 1/2. 

4.3 Solutions using the subalgebra L2 3 

(i) Eq. (J3~7bj) with C ^ 0. If we set A = A p" 2/3 and r] = r] p 2/3 Eq. (j3~7bj) becomes 



2W 



w 2 + 9 -^w 2 + 

2^o 



9a 2 - 1 



-v 



-2 



9s 



2 



W 



-1 



2 ' 2^ 

If we choose r/ = — (|) 1//3 e and a 2 = 1/9, this equation transforms to an equation 
quite similar to PXXXIV 

1 



W 



2W 



W 2 + 4 7 iy 2 - 7]W + 26 2 W~ 1 



where 7 = 3 2 / 3 2- 7 / 3 a A , 5 2 = 3 2 / 3 2- 4 / 3 (C /A ) 2 . We put [EE] 



2 7 W = y + V 2 + - 
2 



and see that V solves 



V = 2V 3 + 7]V + k, 



-- ±4a 5i, 



which is the second Painleve transcendent. Hence we get V = Pji [— (l) 1 / 3 ^ 2 / 3 ] . Since 
W is complex-valued, Ao is a complex constant which must be chosen appropriately to 
make XqW real. 

(ii) Eq. (jSZZEl) with C = 0. We have 



W 



2W 



w 2 + 9 -^w 2 + 



9a 2 -1 



, 9e 



With the choice of a 2 = |, r/ 
PXX 



2r/ 2 

'|) 2 / 3 e and Aq = (I) 5 ^ 3 — , this equation transforms to 



J¥ 2 + AW 2 + 2r]W. 



2W 



Setting U 2 = W we again reach to the Painleve transcendent Pif. 



Hence we write the solution as 

1/2-, 



U = 2U 3 + rjU. 



i> = A ^Pn(7/oP 2/3 ) exp[i(±^ + e lnt + %)}. 



(4.20) 
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4.4 Solutions using the subalgebra L2 4 
Eq. flSJ) for (l3~9b"J) is 

w = ^w 2 -^ + ^ + ±)w + ^- 2 w 2 + ^- 2 w- 

l b a 2 _ 1 A A 2 _ A 
+ f 4( + 2^ _ 32 ~ A£ + 2A 2 ~~ A ' 



(4.21] 



As was done for the subalgebra L 2t i, this equation may be comparable to PXXXIII or 
PXXXIV of Ref. P]. If we compare the coefficients of W 2 and W in (JOT} with (0~7|) 
and ()4.8|) we see that A and 77 must be chosen as 

A = A r 1/3 , r / = r /0 e 1/3 . 

For this choice of A and 77, a comparison of the coefficients with those of the possible 
canonical equations shows us that it can not have the Painleve property. 

If we consider Eq. f!3.9b|) with Co = 0, the same arguments above follows. We stress 
that Eq. (j3.9b|) admits no point symmetries and no constant solutions. 



4.5 Solutions using the subalgebra L2 5 
4.5.1 Equation (I3.11bl) . with C ^ 

By a suitable choice of A and 77 this equation can be transformed to a standard form which 
has no first derivative and a constant coefficient of the quadratic term. The transformation 
and the corresponding transformed equation are given by 

A = A , r) = r) y, 

W = ±-W 2 + i^)W 2 + 2(^) 2 W-\ 
2W Vo A r/ 

In this case we set r? 2 = aoAo/2 and integrate once to obtain the differential equation for 
Weierstrass function or Jacobi elliptic functions (13] 

2C 2 

W 2 = 4(W 3 + CW -5 2 ), C = const, 5 2 = — |. (4.22) 

By virtue of (|4.2jl we require XW to be positive. If we factor the cubic polynomial on the 
right side 

P{W) = A(W 3 + CW - 5 2 ) = 4(W - W X ){W - W 2 )(W - W 3 ), 

where H^ijW^jWa are roots of the cubic polynomial P(W), we see that there are four 
possibilities for the solutions of (|4.22|) : 

(i) P(W) has only one root. This is not possible since P(W) does not include the 
quadratic term W 2 . 

(ii) P(W) has a double root, say W\ = W 2 ^ W 3 and write P(W) as 

P(W) = 4(W 3 + CW - 5 2 ) = A(W - WxfiW - W 3 ). 
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This case is possible for C = —3.2 3 £3 and W^ 2 ,3 are given as W 1>2 = — 2 3^3, W 3 = 
(25)3. Integrating for W > W3 we arrive at the solution given by 

2 1 

M = Cl (-)5(-- + sec 2 r)5, T = Cl y + c 2 , c\ = 3.2- 4 / 3 a A <5 2/3 , 

an o , 

(4.23) 



$ = tan ^y^tanr) - \j^r + $ , *0 = Mexp(i$). 



(iii) P(I¥) has three distinct roots. In this case, roots ^1,2,3 are determined by the 
system 

W 1 + W 2 + W 3 = 0, W X W 2 + W X W 3 + W 2 W 3 = C, W X W 2 W 3 = 5 2 . 

We see that two of the roots have the same sign, while the third is always positive. There 
exist two orderings of the roots yielding different solutions. 

(a) Wi < W 2 < < W 3 (b) < W x < W 2 < W 3 , W% < —. 

Integration for W > W 3 gives both for (a) and (b) 

W = W 2 + (W 3 -W 2 )[cn(r,k)]- 2 , r = c lV + c 2 , k 2 = ^ _ , 

M=(W) 1/2 , c 2 = ^(^ 3 -^i), (4.24) 

5 f dr . , 

$ = ; / — + $ , ib = Mexp(i$) 

with ao, Ao > 0. 

Integrating on W\ < W < W 2 we get 

W = W 1 + (W 2 -W 1 )sn 2 (r,k), r = Cl y + c 2 , k 



, W 2 - W x 



, w 2 - w x 



W 3 -Wi 



M=(W) 1/2 , c\ = a -^{W 3 -W^ (4.25) 

$ = - i [ % + $ , ib = Mexp(i$) 

VW 3 -Wj W y PV ; 

with a , A < for the case (a) and a , A > for the case (b). 

(iv) Finally, we investigate the case where P(W) has two complex roots and a real 
root. Let W 2 = p + iq, W 3 = p — iq, then p, q and W\ are found from the system of 
equations 

2p + W 1 = 0, l- + 2 P W l = C, p 2 + g 2 = l-. 
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Since W\ > 0, integration on W > W\ gives the following solution for do, Ao > 

ur A + W 1 + (W 1 -A)cn(T,k) 

W = l + enfofc) ' t = W + °» 

M = (W) 1 / 2 , e= A + P~W l 



$ = — = [ % + $ , c\ = 2Aa \ 



2A (4.26) 



?/> = Mexp(z$), A 2 = (p- Wx) 2 + q 2 . 

4.5.2 Equation (|3.11hj) . with C = 
From (|4.22j) we have 

W 2 = AW(C + W 2 ), C = const. (4.27) 
Depending on the integration constant C, we obtain the following elliptic solutions for 

(i) C = 0. In this case the solution is written for ao > : 

i> = (~) 1/2 (y - yo)' 1 exp(z$ ). (4.28) 
a 

(ii) C = -p 2 < 0. We write Eq. (PT2TI) as 

W 2 = 4W(W - p)(W + p). 
Integration for W > p and — p < W < gives respectively 



(4.29) 



^ = ci(— ) 1/2 [cn(ciy + c 2 , fc)] l exp(z$ ), a , A > ; 
a 

■0 = c i( — ) 1/2 cn ( c i2/ + c 2, k) exp(i$ ), ao, A < 

a 

with c 2 = a A (-C) 1/2 , A; 2 = 1/2. 

(iii) C = p 2 > 0. The solution for ao, Ao > is given by integrating ()4.27j) : 

2 

V> = ci(— ) 1/2 tn(cij/ + c 2 , fc) dn(ciy + c 2 , fc) exp(z$ ) (4.30) 
a 

with c 2 = aoAl f 1/2 , fc 2 = 1/2. 

4.6 Solutions using the subalgebra L<i$ 
4.6.1 Equation (|3.13hjl . with C ^ 

By a suitable choice of A and rj this equation can be transformed to a standard form which 
has no first derivative and a constant coefficient of the quadratic term. The transformation 
and the corresponding transformed equation are given by 

A = A , r] = rjoy, 
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w = -i-w 2 + (^)w 2 + 2 4w + 2{^-fW~\ 

In this case we set rj 2 . = a X Q /2 and integrate once to obtain 

W 2 = 4(W 3 + aW 2 + CW -5 2 ), C = const., (4.31) 

where 5 2 = 2Cq/ (a \l) and a = 2e/(a \ ). Solution of this equation can be expressed 
in terms of Weierstrass function or Jacobi elliptic functions. By virtue of (J4.2j) we 
require \W to be positive. Let us write the right side as 

P(W) = A(W 3 + aW 2 + CW - 5 2 ) = 4(W - W^iW - W 2 ){W - W 3 ), 

where Wi, W 2 , W 3 are roots of the cubic polynomial P(W). We then have to distinguish 
between four different cases: 

(i) P(W) has only one root. This is the case when 

Wi = ^= (f) 1/2 = 5 2 '\ 



The solution is 

M = (|1 + ^f' 2 , t 2 = Cl y + c 2 , c\ = 8a X , 
Sao r 

35 9 95 275, ,„ . , , , IN , 

$ = — r 2 - — r + — ln(3 + ar) + $ , ip = Mexp[i(et + $) . 
4a 2cr 2cr 

(ii) P(W) has a double root, say W% = W 2 ^ W 3 and write P(W) as 

P(W) = 4(W 3 + aW 2 + CW - 5 2 ) = 4(W - W 1 f{W - W 3 ). 
We find W± and W 3 from the system 

-a = W 3 + 2W X C = 2W 1 W 3 + Wl 5 2 = W 3 W?. 
The roots are ordered as follows: 

(a) Wi > 5 2/3 => < W 3 < W x , (b) W 1 < 5 2/3 => W 1 < W 3 . 
For the case (a) we integrate on W > W\ and W 3 < W < Wl> then find 

W = W x + (Wi - W 3 )cosechV and W = W x + (W 3 - W^sechV 
respectively. The corresponding solutions are 

M=(W) 1/2 , T = c iy + c 2 , c? = ^°(Wi-W 3 ), 

$ = ^rw 3 + % * = Mex ^ £t + $)] ' 



(4.32) 



(4.33) 
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For the case (b), we integrate on W > W 3 and obtain the solution 
W = W 1 + (W 3 -W x ) sec 2 r, T = c lV + c 2 , 

M=(W) 1/2 , c 2 = ^(^ 3 -^i), (4 . 34) 

$ = vm 6 _ m J % + % i> = MeMKet + *)]. 

(iii) P(J¥) has three distinct roots. In this case, roots ^1,2,3 are determined by the 
system 

W l + W 2 + W 3 = -a, W X W 2 + W X W 3 + W 2 W 3 = C, W 1 W 2 W 3 = 5 2 . 

We see that two of the roots have the same sign, while the third is always positive. There 
are two cases to consider: 

S 2 

(a) W x < W 2 < < W 3 (b) < W 1 < W 2 < W 3 , W 2 < —. 

W\ 

Integration for W > W 3 gives both for (a) and (b) 

W = W 2 + (W 3 -W 2 )[cn(r,k)}- 2 , r = c lV + c 2 , k 



,2 w 2 -w l 



W 3 -Wi 

M=(W) 1/2 , C l = ^ { w 3 - Wl ), (4.35) 

5 f dr 



f dr 

= J — + $0, V = Mexp[i(st + $)] 



vW 3 - W^i 

with a , A > 0. Integrating on W\ < W < W 2 we get 

W = W 1 + (W 2 -W 1 )sn 2 (r,k), r = ay + c 2 , k 



, W 2 -W 1 



W 3 -Wi 



M=(W) 1/2 , c 2 = ^(^ 3 -^i), 

with a , Ao < for the case (a) and ao, Ao > for the case (b). 

(iv) Last we investigate the case where P(W) has two complex roots and a real root. 
Let W 2 = p + iq, W 3 = p — iq, then p, q and W\ are found from the system of equations 

2p + W 1 = a, ^- + 2pW 1 = C, p 2 + g 2 = i!-. 

Since W 1 > 0, integration on W > W\ gives the following solution for a , A > 0. 
A + W 1 + (W 1 -A)cn(r,k) 

w = TT^J) ' - = C ^ + C2 ' 

M = (W) 1 / 2 , e= A + p- Wl 



2A ' (4.36) 
Ad 

i> = Mex.p[i(et + $)], A 2 = (p - W x ) 2 + q 2 ■ 



5 f dr 9 

$ = 271 y ^ + $ °' c? = 2Aa ° A °' 
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4.6.2 Equation (13.1 3b jl with C = 

In this case, Eq. for W reduces to 

W = ^-W 2 + (^)W 2 + %W. (4.37) 

We are going to investigate the solutions of Eq. ()4.37|) for e = 1 and e = —1 separately. 
By the choice rj 2 = Ss y 1 = 1 this equation transforms to 

W = — W 2 + AW 2 + 2eW 
2W 

which is PXIX of jTJ for e = 1. The first integral is again the elliptic function equation 

W 2 = 4W(C + eW + W 2 ), C = const. (4.38) 

If we label the right-hand side of Eq. ()4.38|) as P(W) we find the roots of the cubic 
polynomial P(W) to be 



-l-v^4C -l + v / T^4C f 

Wi = , W 2 = 0, Wj, = ^ for e = lj 

and 

^ = 1^4^, = »i = l±4E^ for e = _ L 

We see that ordering of the roots depends on C: 

C < => W 1A < < Wa.e for e = 1, 

1 ' ' (4.39) 

< C < - < W 1A < W 3fi for e = -1. 

Depending on the integration constant C, we obtain the following solutions for 

(i) C = 0. In this case the solution is given below with k = ±1 for suitable sign of 

2 

< V = (— ) 1/2 cosech(/q/ + ci) exp[i(t + $ )] ; 

ao 

-1<W<0, V = ( — ) 1/2 sech(/q/ + ci)exp[i(t + $ )] ; (4.40) 

a 
2 

1 < W, V = (— ) 1/2 sec(/«y + ci) exp[i(-i + $ )] • 

ao 



(ii) C = 1/4. 



< W, ^ = (— ) 1/2 tan(4=2/ + ci) exp[i(t + $ )] 
ao a/2 



0<W<^, i>= (-) 1/2 tanh(^ + Cl )exp[2(-t + $ )] ; (4.41) 
z ao v 2 

i < W, V = (-) 1/2 coth(^y + ci) exp[i(-t + %)] 
l a \/2 
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with a > 0, K = ±1. 

(iii) C < 0. For both e = ±1, the roots of P(W) have opposite signs. An integration 
with respect to W 3i6 < W and < W < W 2 $ results in following solutions: 

V> = ( V/1 ~ 4C ~ £ ) 1/2 [cn(ci|/ + c 2 , ft)]" 1 exp[i(et + %)], a > ; 

gjL + (4.42) 

V> = ( -) 1/2 cn(dy + c 2 , k) exp[i(et + $ )], a < 

— a 

with c 2 = VT^4C, fc 2 = ^gf. 

(iv) < C < j. We list the solutions following the interval that the integration has 
been taken on: 

W 2 < W, 

^ = ci( — ) (- = tn(c 1 y + c 2 ,k)exp[i(t + ® )}, a > ; 

Va / V1 + V1-4C/ f443 , 

w x < w < w 3 , y ' 1 

i> = ci(— ) l/2 [ ]~ lZ^ + cn 2 { Cl y + c 2 ,k)} 1/2 exp[i{t + %)}, a < 
\ Qo 7 1 + yl ~ 4G 

with r 2 - ^v 71 ^ b 2 - 

witn q - 2 , ft - 1+v/I= 4c- 

W 6 < W, 

1> = ( 1 ~ V/I ^ ) 1/2 [1 + -^^^(cn(c l2/ + c 2 , A:))- 2 ] 1 / 2 eX p[i(-t + $ )] ; 

a 1 - VI -4C (4 44) 

W 5 < W < W 4 , 

,j, = ( Iz 1 ~ AC yi 2 sn(c iy + ca, fc) exp[i(-t + $o)] 
ao 

witha >0, c? = fc2 = Wgg. 

(v) T < C. We have two solutions for a > : 



^ = (— ) 1 / 2 C 1/4 tn(C 1/4 i/ + c h ft) dn(C 1/4 y + c x , ft) expfzfrf + $„)] (4.45) 
a 



with ft 2 = 



4.7 Solutions using the subalgebra L2 7 

This subalgebra leads to second order non-linear ODEs of which solutions are written in 
terms of Painleve transcendents. 

(i) Eq. (j3.15b|) with Cq 7^ 0. If we set A = Ao, r\ = ??o£, this equation transforms to 
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For the choice r] = — 1 this equation becomes 

W = t^ttW 2 + -rjW + 25 2 W-\ 

— V V 

where 7 = aoAo/2 and S 2 = (C /A ) 2 . This equation is quite similar to equation PXXXIV 
of jTT] and is integrated in the same way. We put 

2^W = V + V 2 + - 
2 

and see that V is solved by the equation 

V = 2V 3 + rjV + k, k = -~±4a 5i, 

which is the second Painleve transcendent. Hence we find V = Pn(—£)- Since W is 
complex-valued, Ao is a complex constant which must be chosen appropriately to make 
X W real. 

(ii) Eq. (j3.15b|) with Co = 0. Under the same transformation X = X , i] = rjoC, this 
equation reduces to 

w = — w 2 + ( 2 -^)w 2 + ±w. 

2W 1 r, 2 } V 3 

With choice r] = 2~ x l 3 and Xq = 2 1 / 3 /a it is nothing but the canonical Eq. PXX of JT] 



W = ^W 2 + AW 2 + 2r]W. 

If we set U 2 = W, the second Painleve transcendent Pjj makes its appearance again 

U = 2U 3 + 7]U. 

We can write the solution as 

V> = ^Pii(2^ 1/3 exp[i(| -£ + $ )]. (4.46) 

4.8 Solutions using the subalgebra L210 
Eq. (gH) for (EHTEl) is 

w = ±-w 2 - i(? + ±)w + 2 -^w 2 + 2{% 2 W~' 

2W rj i] X i] A Xr\ 



(4.47) 



Just as for the subalgebra L 2i i, this equation seems a candidate to be transformable to 
PXXXIII or PXXXIV. By comparison of the coefficients of W 2 and W in 1EQ7I) with 
(|4.7|) and ()4.8|) it turns out that A and 77 must be chosen to be 

A = A , 77 = ??o£- 
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For this choice of A and rj ()4.47|) becomes 



W 



2W 



W 2 + ^0 W 2 + 21 Co v2 



v'o 



W- 1 + (\r? - ^) W. (4.48) 



Since the coefficient of W includes the term i] 2 , Eq. ()3.21b|) cannot be transformed to one 
of the canonical forms, thus it does not have the Painleve property. 

If we consider Eq. (J3.21bj) with Co = 0, the same arguments above apply. We finally 
state that Eq. (J3.21bj) admits no point symmetries and no constant solutions. 



4.9 Solutions using the subalgebra L2, 
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Singularity analysis and integration of Eq. (|3.25b|) are performed in Ref.jH] in which the 
authors showed that Eq. ()3.25b|) for C = 0, namely 

2*X - X 2 + ( j - 2a)x 2 - \x 2 ~ 4a X 3 = 
is of Painleve type and allows a first integral of the form 



V 



■ ^(X - XO 2 + 2«X 2 + 2a + X 3 + C X - 



They also showed that the solution of (j4.49|) is given by 

'w 2 



x(0 



k 

a 



AW 



^_^ + ( W + a) ^ + A + 2 (a + 1) , + _ 



where W(z) satisfies the Painleve transcendent Pjy 
W 



^W 2 + 3 -W* + 4 Z W 2 + 2( Z 2 -a)W + ^, 



where z = /c^, k 4 = —1/16 and the constants a, f3 are determined by 

-2(a + l) 2 + 3/? = 4(2a 2 - 3Ca ), 

(a + l)[2(a + l) 2 + 9(3} = (a/A; 2 ) (4a 2 - 9Ca ). 

To sum up, the solution to (jl.lj) is given by the formula ()3.24|) with 

£(1 + *) 



M 



X- 



4x 



-de 



(4.49) 



(4.50) 



(4.51) 



Remark: We have not been able to successfully treat the reduced equation corresponding 
to the subalgebra L211. 



4.10 Invariant solutions from three-dimensional subalgebras 

We conclude by looking for further solutions on the basis of three dimensional subalgebras. 
Invariance under those subalgebras gives rise to algebraic equations only rather than 
ODEs. We simply list nontrivial solutions in table EU ruling out those leading to trivial 
ones. 
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5 Concluding remarks 



In this paper we considered a cubic Schrodinger equation in 2+1- dimension and gave 
a comprehensive investigation of group-invariant solutions by performing reductions to 
ODEs using appropriate subgroups of the symmetry group and by solving them when 
possible. Our approach is a powerful combination of the group theory and singularity 
analysis. In this respect, it is very similar to the one used in a series of papers for 
an investigation of group-invariant solutions of a 3+1-dimensional generalized nonlin- 
ear Schrodinger equation (which, as special cases, contains dissipative, cubic and quintic 
Schrodinger equations) in PQ 121 El- The solutions obtained represent a conjugacy class 
of solutions. The entire class can be produced by applying a general symmetry group 
transformation ()2.10|) - (j2.11[) to the representative solution. The solutions include Jacobi 
elliptic functions, their degenerate forms as elementary ones and Painleve irreducibles. 
What emerges from a comparison of results in 3+1- and 2+1- dimensions is that there are 
reductions to ODEs with similar structure in both dimensions (for instance, Painleve tran- 
scendents are ubiquitous) while the corresponding invariant solutions have truly different 
expressions. 

We mention that some similarity solutions of (jl.l|) have already been discussed in 
[Tl] in which the author merely chooses an arbitrary combination of the generators for 
performing a symmetry reduction. Reduction to an ODE is achieved in two successive 
steps, first to a 1+1-dimensional PDE and then to an ODE by imposing an inherited 
symmetry of the reduced PDE. On the contrary, in the present article we classify solutions 
into a list of conjugacy classes by using a rigorous subgroup classification so that every 
solution is conjugate to precisely one in the list. 
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